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I .  INTRODUCTION 


Because  of  their  very  large  volume  and  extreme  sensitivity,  extrinsic 
infrared  detectors[l]  operating  above  the  atmosphere  may  be  quite  vulnerable 
to  single  ionization  events  resulting  from  energetic  ionized  particles  in 
space.  This  paper  describes  a  method  of  applying  existing  techniques  for 
calculating  the  intensity  distribution  of  false  signals  resulting  from  cosmic 
rays  as  they  impinge  upon  a  silicon  detector  of  given  dimensions.  The  result 
is  a  plot  of  the  amplitude  distribution  expected  from  near  earth  cosmic  ray 
events.  The  program  can  also  make  comparisons  of  various  shielding  models  and 
the  effects  of  different  particles  and  particle  spectra. 


II.  SPHERICALLY  SYMMETRIC  SHIELDING  MODEL 


Figure  1  Illustrates  how  the  "straight  ahead"  approximation  is  used  to 
determine  the  particle  flux  inside  a  spherical  cavity.  Data  from  range-energy 
tables  are  functionalized  by  means  of  an  interpolation  program  yielding  R(E), 
range  as  a  function  of  energy,  and  its  inverse,  E(R).  Then 

Ej  =*  E(R(E2)  +  Ts)  (1) 

where  Ej  is  the  energy  of  the  particle  incident  on  the  shield  and  E2  is  the 
energy  of  this  particle  after  it  is  degraded  by  the  shield  thickness  Tg. 
Equation  (I)  is  also  the  basis  for  determining  deposited  energy  in  the  sensi¬ 
tive  region  of  the  detector.  The  degraded  flux  spectrum,  which  for  spherical 
geometry  is  assumed  to  be  the  flux  incident  on  the  detector,  becomes 

W  =  *1  (V  x  <dVdV  (2) 

where  <D^  (E^)  is  the  flux  spectrum  incident  on  the  shield  as  obtained  in  this 
case  from  the  work  of  Adams  et.  al.(2],  and  the  derivative  (dE^/dE2)  is  ob¬ 
tained  from  Eq.  (1). 

The  effects  of  shielding  on  the  proton  and  helium  spectra  are  illustrated 
in  Figures  2  and  3,  respectively,  for  two  widely  different  shielding  thick¬ 
nesses.  Note  that  once  the  very  low  energy  component  is  removed,  the  shield 
thickness  has  little  effect  on  the  spectrum. 

Secondaries  coming  from  the  shield  were  neglected  in  this  calculation. 
There  are  two  sources  of  secondaries  to  be  considered:  (I)  those  resulting 
from  nuclear  reactions  between  the  incident  energetic  nuclei  and  the  nuclei  of 
the  shield  material,  and  (2)  those  resulting  from  hot  electrons  produced  by 
the  highly  ionized  tracks  of  the  incident  nuclei.  The  cross  sections  for 
nuclear  reactions  are  known  and  indicate  this  source  of  secondary  generation 
to  be  negligible  compared  to  the  primary  particle  event  rate.  No  attempt  was 
made  to  estimate  the  secondaries  resulting  from  hot  electrons  but  it  seems 
reasonable  that  such  events  if  present  would  produce  only  low  energy  or  low 
magnitude  pulses. 
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WHERE:  Ej  •  E  ( R I E^)  +  T$) 


Figure  1.  Spherically  Symmetric  Shielding  Model 


III.  DEPOSITED  ENERGY 


The  omnidirectional  particle  flux  inside  the  shield  creates  tracks  of 
highly  ionized  regions  (hole-electron  pairs)  in  the  detector.  Essentially  all 
of  the  energy  deposited  within  the  sensitive  volume  of  the  detector  creates 
hole-electron  pairs  (one  pair  per  3.6  eV)  which  are  collected  and  amplified  by 
the  electronics  to  produce  an  output  pulse.  To  get  a  pulse  height  distribu¬ 
tion  we  first  determine  the  deposited  energy  distribution. 

Consider  the  minimal  event  which  we  define  as  the  passage  of  a  proton  of 
energy  ELin,  corresponding  to  the  minimum  in  (dE/dx) ,  perpendicularly  through 
the  smallest  dimension  of  the  parallelepiped  shaped  sensitive  volume.  Then 
the  resulting  deposited  energy  corresponds  to  the  energy  of  the  maximum  of  the 

O 

deposited  energy  histogram.  For  protons  dE/dx  at  Emin  is  1.7  MeV-cm4/gm.  The 

deposited  charge  in  Coulombs  is  (dE/dx  x  10^  x  s  x  1.6  x  10~^/3.6)  = 

(7.56  x  10‘14  x  s)  where,  in  this  case,  s,  the  chord  length  in  gm/cm^,  is  the 

thickness  of  the  detector.  For  a  100  um  thick  (0.0233  gm/cm^)  extrinsic 

detector,  the  minimal  event  deposits  a  charge  of  1.76  x  10~^  pC,  or  about  10^ 
e.h.  pairs.  This  is  quite  large  compared  to  many  infrared  signals  likely  to 
be  encountered.  Assume  the  pulse  is  amplified  by  a  conventional  trans¬ 
impedance  amplifier  (TIA)  as  in  Figure  4.  Since  the  pulse  is  very  fast,  it 
can  be  shown  that  the  output  of  the  TIA  is  VQ  =  Q/Cpg,  i.e.,  the  deposited 
charge  divided  by  the  feedback  capacitance.  For  C^g  *  .01  pf  VQ  is  about  0.2 
volt,  much  larger  than  the  usual  sources  of  noise. 

Calculations  of  single  event  rates  for  microelectronics  invariably  assume 
that  the  deposited  energy  is  (LET)  x  (chord  length,  s,  of  the  track  through 
the  sensitive  volume).  LET,  the  linear  energy  transfer,  is  dE/dX.  This  "LET 
approximation"  is  quite  safe  for  microelectronics  because  of  the  very  small 
dimensions  involved.  For  IR  detectors,  in  some  cases  hundreds  of  microns 
thick,  the  LET  approximation  may  well  be  inadequate.  Thus,  for  the  present 
calculations,  the  average  deposited  energy  DE  is  computed  from  the  range- 
energy  relation  for  each  component  of  the  cosmic  background: 


Figure  4.  Typical  Detector  with  Trans-Impedance  Amplifier 


E2, 


for  E2  <  E(s) 


DE 


E2  -  E(R(E2)-s),  for  E2  >  E(s) 


(3) 


where  E2  is  the  energy  of  the  particle  after  it  has  been  degraded  by  the 
shield,  but  before  it  has  entered  the  detector.  Note  that  the  total  incident 
energy,  E2,  is  deposited  in  the  detector  if  it  is  less  than  the  energy  re¬ 
quired  for  a  range  of  s.  If  E2  is  greater  than  this  the  particle  will  escape 
the  sensitive  region  and  DE  is  E2  minus  the  energy  of  the  particle  as  it 
escapes  from  the  sensitive  region.  To  account  for  collection  of  charge  by 
diffusion  one  may  add  a  small  amount  to  the  dimensions  of  the  detector  to 
estimate  its  effective  dimensions.  This  has  not  been  done  in  the  present 
calculations. 

Figures  5  and  6  plot  DE  (E2)  for  two  widely  differing  values  of  s:  0.1 

•y 

and  1.75  gm/cm  ,  respectively.  Also  plotted  with  the  dashed  line  is  the  LET 
approximation  (s  x  dE/dx).  It  is  seen  that  the  LET  approximation  is  very  good 
for  those  protons  that  are  not  stopped  within  the  sensitive  volume.  But  for 
those  stopped  (E2  <  E(s))  the  LET  approximation  greatly  overestimates  the 
deposited  energy.  For  s  =*  0.1  gra/cra  ,  there  are  very  few  protons  of  energy 
E2  <  E(s)  =  6.6  Mev  (.02%).  For  s  *  1.75  gm/cm^  about  0.2%  of  the  incident 
protons  are  less  energetic  than  E(s)  ■  4.0  MeV  (see  Figure  6).  Thus,  using 
the  LET  approximation  could  affect  the  shape  of  the  pulse  height  histogram  in 
the  very  large  pulse  height  region  if  the  dimensions  of  the  sensitive  volume 
are  large  enough.  If  the  maximum  chord  length  is  less  than  about  0.75  cm 
(1.75  gm/cm4)  the  LET  approximation  should  be  reasonably  accurate.  In  the 
present  calculations,  however,  the  more  accurate  method  (Eq.  (3))  is  retained 
for  purposes  of  generality. 


IV.  AMPLITUDE  DISTRIBUTION  FUNCTIONS 


The  function  DE  (E2)  for  a  given  chord  length  s  through  the  sensitive 
volume  can  always  be  divided  into  three  monotonic  intervals:  (0  <  E2 
<  E(s) ) ,  (E(s)  <  E2  <  ^rain)*  and  (E2  >  ^inin^  *  Figures  5  and  6  the  minimum 

in  DE  corresponds  to  Emin.)  Let  Ng  (DE)  be  the  deposited  energy  distribution 
for  a  given  s;  i.e.,  Ng  (DE)  *  ADE  is  the  number  of  events/cm^-ster  that 
deposit  energy  in  the  range  DE  to  (DE  +  ADE).  Thus,  the  contribution  to 
Ng(DE)  by  the  ith  interval  is: 


N  (DE)  x|ADE| 


2U 

/  $?(E  )  dE, 

E„. 


where  EoT  and  E?I,  are,  respectively,  solutions  of  Eq.  (3),  as  follows: 


DE  =  DE  (E2L,s) 


DE  +  ADE  »  DE  (E,^,  s)  (5) 

Equation  (4)  assumes  that  the  probability,  P(DE),  of  depositing  energy  DE  is  a 
delta  function:  P(DE)  ■  6  (DE  -  DE).  This  is  a  very  safe  assumption  for 
protons  and  heavier  particles.  The  total  contribution  is  the  sum  of  the 
three;  N„(DE)  =  EN  .(DE)  provided  the  ADE  intervals  are  the  same.  For  the 
total  distribution  Nj(DE)  of  deposited  energy  produced  by  the  jth  component  we 
utilize  the  exact  chord  length  distribution  of  Petroff[3J,  Fg,  as  follows: 


N  (DE)  »  4ir  x  A  x  /  N  (DE)  x  F  ds 


where  A^  is  the  average  projected  area  of  the  sensitive  volume  and  Smax  is  its 
largest  chord  length: 


(L2  +  W2  +  H2)0,5 


where  L,  W,  H  are,  respectively  the  length,  width  and  height  of  the  parallele¬ 
piped  shaped  sensitive  volume.  The  total  distribution  is  the  sum  of  those 
produced  by  individual  components:  N(DE)  =  J  N^DE).  In  the  present  calcula¬ 
tion  we  considered  protons  and  helium  only.  In  other  calculations  we  found 
that  the  inclusion  of  heavier  elements  make  no  noticeable  difference.  The 
cumulative  distribution  is 


CN(DE) 


DE<Smax> 


/  N(DE)  x  dDE 
DE 


(8) 


One  can  calculate  the  total  number  of  cosmic  ray  particles  per  cm  of 
average  projected  area  that  penetrate  the  sensitive  region  by  integrating 
Adam's  spectra  (after  a  moderate  amount  of  shielding).  This  gives  1.64 
events/cm^sec  and  agrees  quite  well  with  CN(O)  obtained  from  Equation  (8), 
confirming  the  internal  consistency  of  the  calculation. 

To  convert  Eq.  (8)  to  output  voltage  distribution,  N(V)  =■  N(DE)  x  ) 

where  the  peak  output  voltage,  V,  is  obtained  from  linear  circuit  theory  for 
the  circuit  of  Figure  4: 


V 


V 

CFB  norm 


(9) 


where  Cpg  is  the  capacitance  across  the  feedback  resistor,  RpB,  of  the  TIA,  G 

is  the  D.C.  gain  of  the  following  amplifiers,  and  Q  is  the  collected  charge: 

Q  »  4.44  x  10“^®  x  (DE).  From  linear  circuit  theory,  Vnorm  can  be  determined 

as  a  function  of  xD,  x_,„,  and  xA,  which  are  respectively  the  response  time  of 

the  detector  (carrier  lifetime),  RpB  CpB,  and  I / 2xxf q ,  where  fQ  is  the  roll  off 

frequency  of  the  following  amplifiers.  If  x^  is  very  small  compared  to  the 

smaller  of  x  and  x.  (which  is  usually  the  case)  then 
r  d  A 

_  _-[ z/ ( z— 1 ) ] 


norm 


(10) 


where,  Z  *  The  derivative  — rr  Is  obtained  from  Eq.  (9). 

A  rB  DE; 

Equation  (10)  is  plotted  in  Figure  7  to  show  its  monotonic  behavior. 

Thus  to  decrease  the  effect  of  the  spurious  cosmic  ray  pulses  one  should 
increase  the  amplifier's  response  time  and  decrease  the  parasitic  capacitance 
across  the  feedback  resistor  of  the  TIA  as  much  as  possible  without  loss  of 
signal  amplitude.  This  requires  keeping  less  than  the  dwell  time  or  the 
pulse  width  of  the  optical  signal. 

In  Figure  8  the  cumulative  amplitude  distribution  calculated  for  an  array 
of  extrinsic  detectors  shielded  with  62  gm/cm^  of  high  Z  material  is  plotted 
as  a  function  of  amplitude  in  arbitrary  units.  (An  array  of  N  identical 
detectors  would  have  an  event  rate  of  N  times  that  of  a  single  detector.) 

This  array  has  also  been  flight  tested  at  low  orbit  where  Van  Allen  radiation 
can  be  neglected,  but  above  the  earth's  atmosphere.  Thus,  a  likely  candidate 
for  the  cause  of  fast,  spurious  output  pulses  would  be  near  earth  cosmic 
rays.  The  effect  of  earth  shielding  was  taken  into  consideration  by  assuming 
exposure  to  only  2x  steradians  of  omnidirectional  cosmic  rays.  The  cumulative 
amplitude  distribution  of  "fast,  spurious"  output  pulses  obtained  during  the 
flight  test  is  also  plotted  in  Figure  8  against  the  same  scale.  Considering 
the  uncertainty  in  particle  fluxes,  the  agreement  is  very  good,  except  at  low 
pulse  heights  where  the  calculated  values  are  too  low.  This  might  be  the 
result  of  secondaries  produced  in  the  thick  shield  since  these  secondaries 
were  not  included  in  the  calculation  and  are  expected  to  lie  in  the  low  ampli¬ 
tude  region. 


V.  CONCLUSIONS 

Good  agreement  was  obtained  between  the  calculated  amplitude  distribution 
of  cosmic  ray  events  and  the  spurious  events  observed  in  an  extrinsic  detector 
during  a  flight  test  below  the  Van  Allen  radiation  belts.  The  date  of  the 
flight  was  taken  into  consideration  according  to  the  recipe  described  by 
Adams [2].  Although  the  deposited  energy  was  determined  from  the  range-energy 
relation  it  appears  that  the  LET  approximation  would  be  adequate  for  detectors 
whose  maximum  dimension  is  less  than  0.75  cm.  A  method  of  decreasing  the 
amplitude  of  these  pulses  by  taking  advantage  of  their  fast  risetime  was 
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Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy -i nduced 
env i ronments . 

Space  Sciences  Laboratory :  Magnetospherlc ,  auroral  and  cosmic  ray 
physics,  wave -part ic le  interact  ions ,  magnetospherlc  plasma  waves;  atmospheric 
jii>i  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


